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Abstract 


Among natural hazards, floods are of special importance in terms of financial and human losses. 
According to available statistics and information, flood damage in most parts of the world, 
especially in developing countries, is the the most distructive and generally makes socio- 
economic and environmental damages. So, efforts for reducing the damage are the duties of 
related agencies. Estimating and determining the maximum amount of flood is one of the first 
steps for reducing its damage. Therefore, providing a definition of it is a fundamental step. In 
this regard, the selection of the maximum possible flood plays an important role. By definition, 
the maximum possible flood is the largest flood that can reasonably be expected to occur in a 
particular region and time (World Meteorological Organization, 1986). Frequency analysis is 
one of the reliable methods in flood estimation if there are sufficient statistics and it is possible 
to fit a suitable distribution function. This method, which uses a complete statistical series (all 
available data), often uses normal statistical distributions, normal logs, three-factor normal logs, 
limit values distribution, Pearson, Pearson logs, and WebStud program (13). The purpose of 
frequency analysis of events in hydrology is to obtain the maximum probability of occurrence of 
values. The selected area for studying the geomorphic basin of Shiraz is in the sections of 
Aliabad, Maroon, Eghbalabad, Chenar Sukhteh and Bagh-e-Safa stations. In this study, long- 
term precipitation statistics in the rain gauge station of Shiraz geomorphic basin over a period of 
50 years and maximum daily discharge statistics in a period of 44 years were used. 
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حغرافیا و مخاطرات محیطی. سال دوازدهم. شمارة جهل و هشتم. زمستان ۰۱۶۰۳ صص ۱۱۷-۱۳۷ 
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حکیده 


پژوهش حاضر با هدف شناسایی مطالعه جامع» جهت مقادیر دبی‌های موثر بر وقوع سیلاب‌های شهری 
رودخانه خشک شیراز انجام شده است. حوضه آبریز شیراز با مساحت برابر ۱۸۱۵/۱۰ کیلومترمربع و اقلیم 
خشک و نیمه‌عشک» همراه با پوشش گیاهی متوسط به دلیل خشکسالی‌ها و تغییرات کاربری اراضی است 
که نشان از شرایط سیلابی حوضه در زمان رخداد رگبارهای ناگهانی دارد. روش تحقیق آماری- تحلیلی و 
نوع ol‏ کاربردی است. بدین منظور در جمع‌آوری داده‌ها از آمار درازمدت بارش در ایستگاه باران‌سنجی 
حوضه آبریز شیراز در یک دوره ۰ ساله و آمار دبی حداکثر روزانه در یک دوره ۶۶ ساله ایستگاه‌های 
آب‌سنجی موجود در منطقه. پس از بررسی موقعیت مکانی» کیفیت تجهیزات. مدت و کیفیت داده‌های 
ایستگاه‌های آب‌سنجی, استفاده شده است. سپس به‌منظور انتخاب توزیع مناسب: داده‌های هر ایستگاه در 
محیط "۳616" وارد شده و اقدام به استخراج ماکزیمم مقادیر دبی شده و خروجی‌ها و ضوابط جبری 
در نرم‌افزار 16 Graphers‏ محاسبه و تهیه شده است. در مرحله بعد جهت رفع اوریب يا بایاس از 
میانگین و انحراف استاندارد مقادیر دبی در تابع احتمال تیپ یک گمبل» با محاسبه گشتاور دوم توزیع حد 
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آستانه از ساعتار جز به جز ریمان استفاده شده است. درنهایت این تحقیق نشان داد تا قبل از قاجاریه Lids‏ 
توسعۀ فیزیکی شهر شیراز خارج از محدوده عبور رودخانه خشک بوده است؛ اما به‌تدریج با روند توسعه 
شهرنشینی» عدم رعایت حریم رودخانه و مدیریت‌های غلط شهری در حال حاضر کوچک‌ترین بارش‌های 
سالانه موجب وقوع سیلاب شهری شیراز شده است؛ به‌طوری‌که با این تفاسیر بروز فاجعه در به وقوع 
پیوستن دبی پیک طی صدسال با احتمال ۷۰ در صد با حدود ۲۰۰ مترمکعب در ثانیه حتمی است. اما این 
به معنای وقوع سیل در شیراز طریق رودخانه خشک نبوده است. 

کلیدواژه‌ها: سیل. دبی حداکثر: توزیع گمبل. حوضه شیراز. 


۱- مقدمه 


سیل ممکن است در اثر بروز یک بارندگی شدید و یا مستمر و یا ذوب ناگهانی برف‌ها در یک حوضه آبخیز به 
وقوع بپیوندد (اصغری مقدم» ۱۳۸۲). پهنه‌بندی سیلاب تقسیم‌بندی سطح زمین به یک سری از پهنه‌های دارای ارزش 
معین است. پهنه‌بندی سیلاب به تعیین ناحیه‌هایی در داخل دشت سیلابی GUE!‏ می‌شود که برای کاربری‌های مختاف 
از قبیل فضاهای باز تفریحیء کشاورزی, محوطه‌های صنعتی و مسکونی مورد استفاده قرار می‌گیرد و همچنین با 
توجه به تعاریف وزارت نیرو پهنه‌بندی سیل را این چنین تعریف می‌کند: تعیین پهنه‌های سیل گیر برای سیلاب‌های با 
دوره بازگشت‌های معین را پهنه‌بندی سیل می گوییم. (نشریه ۳۰۷ وزارت نیرو). مخاطرات طبیعی به‌عنوان تهدیدی در 
حال رشد در ارتباط با رفاه و توسعه جوامع مطرح‌اند و حسارت‌های ناشی از مخاطرات طبیعی» هم از دیدگاه تکرار 
و هم از نظر صدمه‌هایی که به بار می‌آورند» بر جامعه اثر می‌گذارند (پورطاهری و همکاران؛ ۱۳۹۰). تشدید صعودی 
خسارت سیل در دو دهه گذشته سبب شده که آرزوی دیرینه درباره حل قطعی مسئله سیل و رواناب‌ها جای خود را 
به واقع‌گرایی و درک این واقعیت دهد که همیشه نمی‌توان در مهار سیلاب‌ها موفق بود بلکه باید کوشید تا پیامدهای 
مخرب و زیان‌بار آن را کاهش داد (صفاری و همکاران» ۱۳۹۰). در این میان رودخخانه‌ها به‌عنوان اصلی‌ترین منابع 
تأمین کننده آب برای انسان و سایر موجودات به شمار می‌روند و بعضاً این منشاً زندگی باعث نابودی و وارد شدن 
Chel‏ جرا فار می شود کار این تشن رفن مرو یکی روتسا دز مقایل بت‌لات‌های اما 
برای کاهش حسارت‌های وارده بر مناطق شهری و روستایی» تأسیسات در حال ele‏ مزارع و ple‏ کاربری‌های 
موجوه در اطراف رر ههه دازای اهمیت ead) th, po ely‏ و کر CAPE‏ زیرا ر آن‌ها Ble. ge‏ 
تصمیم‌ها و راهکارهایی را برای مهار سیلاب و به حداقل رساندن خسارت‌های ناشی از آن کسب نمود (حجازی و 
همکاران, ۱۳۹۸). در ادامه به بررسی برخی از تحقیقات داخلی و خارجی انجام گرفته در مورد حطر سیلاب پرداخته 
شده است. اهداف پهنه‌بندی با توجه به نیازها می‌تواند متفاوت باشد تعیین بستر و حریم تهیه نقشه خطر سیلاب» 


سال دوازدهم بررسی حداکثر دبی مقادیر سیلاب حوضه آبریز ... na‏ 


ارزیابی خطر سیلاب انجام گرفته است که در ادامه به چند مورد آن اشاره می‌شود: (حسین‌زاده و همکاران, ۱۳۹۱) به 
شبیه‌سازی سیلاب رودخانه کشکان به مطالعه پرداخته‌اند نتایج حاصل از پژوهش OLS‏ می‌دهد که روند افزايش سطح 
سیل گیر در دوره بازگشت‌های ۵ تا ۱۰۰ سال سیر صعودی‌تر دارند؛ درحالی که مساحت پهنه سیلاب در دوره 
بازگشت ۱۰۰ ساله و ۲۰۰ ساله اختلاف چندانی باهم ندارند. در ایران حدود۲ درصد اراضی مستعد ابتلا به مخاطره 
سیلاب هستند و 3 درصد دارای احتمال خطر سیلاب با ریسک بسیار زیاد است (خسروی و همکاران. ۲۰۲۱). در 
سال‌های اخیر ارزیابی‌ها و برآورد سیلاب به دلیل افزایش شهرنشینی و تمرکز سکونتگاه‌های مسکونی در اطراف 
رودخانه‌ها و دشست‌های سیلابی به‌طور قابل‌توجهی افزایش adh‏ است(آوند! و همکاران ۲۰۲۱). سیل یکی از 
مخرب‌ترین مخاطرات طبیعی است؛ بخصوص در مناطق خحشک و نیمه‌عشک L)‏ سیل‌های ناگهانی ) که بارند گی 
زیادی را تجربه می‌کنند. این مخاطره ناشی از پیوند چندین عامل طبیعی و انسانی است. سیلاب هنگامی رخ می‌دهد 
که دبی رودخانه به حدی افزایش می‌یابد که آب از سواحل رودخانه عبور کرده و مناطق مجاور را غرقاب می AS‏ 
درنتیجه جوامعی را که در مجاورت رودخانه‌ها ساکن هستند. تحت تأثیر قرار می‌دهد علاوه بر آن تغییرات آب و 
هوایی احتمال وقوع سیل و بزرگی آن‌ها را افزایش می‌دهد (ترمه "و همکاران» ۲۰۱۸ ؛ ازین "و همکاران ۲۰۲۰). 

سیل سالانه به‌طور متوسط VO‏ میلیون نفر را در سراسر جهان تحت تأثیر قرار می‌دهد (حدود ۲۰۰۰۰ تلفات در 
هرسال)» ضرر اقتصادی حاصل از آن ۹۵ میلیارد دار پوده است در سال‌های اخیر فراوانی سیل بیش از ۸۶۰ افزایش 
یافته است؛ زیرا تغیبرات آب و هوایی به‌طور قابل‌توجهی بر شدت الگو و بزرگی‌های سیلاب‌ه تأثیر داشته است. 
به‌ویژه کشورهای آسیایی با برخورد نسبتاً زیاد در معرض سیل و زلزله قرار دارند در آسیا بیش از ۸٩۰‏ تلفات انسانی 
ناشی از بلایای طبیعی به علت جاری شدن سیل است (جانک* و چن ۲۰۱5؛ خسروی" و همکاران ۲۰۱۸ ۲۰۲۱). 
oly!‏ اخیراً مسیلاب‌های فاجعه باری را تجربه کرده اسست که علاوه بر تأثیر تغییرات آب و هوایی» تأثیر فعالیت‌های 
انسانی ناشی از مدیریت ضعیف حوضه آبریز ازجمله جنگل‌زدایی. چرای بیش‌ازحد دام. توسعه سکونتگاه‌های 
روستایی و شهری باعث کاهش clad‏ هیدرولیک رودخانه‌ها شده است» سیل‌های مارس و آوریل ۲۰۱۹ در ايران که 
۵ استان از ۳۱ استان کشور را تحت تأثیر قرار داده است درنتیجه بیش از ۷۷ نفر تلفات انسانی و خساراتی به ارزش 
۲ میلیارد دلار به سکونتگاه‌های شهری و روستایی و زمین‌های کشاورزی وارد شده ات 
-(https://en. wikipedia.org/wiki/2019_Iran_floods)‏ در پژوهشی دیگر پهنه‌بندی سیلاب با استفاده از مدل 


هیدرولیک CHEC-RAS)‏ در رودخانه سرباز استان سیستان و بلوچستان به مطالعه پرداخته شده و ابتدا دبی با دوره 
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باز گشت‌های مختلف با استفاده از مدل SCS‏ در نرم‌افزار 11۳0-11۷15 محاسبه شده و نقشه پهنه‌بندی سیلاب حوزه 
آبخیز سرباز محدوده شهر راسک با دوره بازگشت‌های مورد نظر با کمک مدل هیدرولیکی HEC-RAS‏ تهیه شده 
است (دامادی و همکاران, ۱۶۰۰). SULE)‏ و همکاران (TVW‏ به تهیه نقشه دشت سیلابی با کمک HEC-RAS‏ و 
GIS‏ رودخانه کابل پاکستان به مطالعه پرداخته‌اند که در این مطالعه با کمک داده‌های هندسی مدل رقومی ارتفاعی 
DEM‏ داده‌های هندسی رودخانه تهیه شده و دوره بازگشت‌ها را بیان کرده است که ۶0 درصد ناحیه در معرض سیل 
قرار دارد که با تصویر MODIS‏ مقایسه شده و نتایج مدل مورد تائید واقع شده است. (کابنگی و همکاران ۲۰۱۷۰) به 
توصیف مخاطره سیلاب در مناطق فاقد اطلاعات با استفاده از تکنیک سنجش ازدور و شاخص احتمال خطر سیلاب 
مبتنی بر GIS‏ به مطالعه پرداخته‌اند؛ فاکتورهایی از قبیل شیب. جریان تجمعی, تراکم شبکه زهکشی فاصله از SES‏ 
زهکشی کاربری اراضی شدت بارندگی و زمین‌شناسی به‌عنوان عوامل ایجادکننده سیل مورد استفاده قرار گرفت که 
همه این فاکتورها از داده‌های سنجش زدور LANDSAT SRTM‏ به جز زمین شسناسی استخراج گردیده است و با 
کمک آنالیز سلسله مراتبی SE‏ فاکتورها مشخص شده و درنهایت نقشه سیلاب تهیه و کاربری متناسب با مناطق سیل 
پیشنهاد شده است. (زلناکووا" و همکاران ۲۰۱۹) به مدل‌سازی احتمال خطر سیلاب رودخانه -اسلاوینس در 
اسلواکی به مطالعه پرداخته‌اند که با تلفیق 618 و ارزیابی آسیب‌های سیل را با استفاده از تحلیل هزینه (CA)‏ انجام 
داده‌اند که این مدل در (HEC-RAS)‏ اجرا شده است که درنهایت مناطق احتمال خطر سیلاب مشخص گردیده است 
و میزان حسارت محاسبه شده است. در تمامی این تحقیقات به مطالعه پهنه‌های سیلابی. دوره‌های بازگشت مختلف و 
polis‏ آب‌گرفتگی در هر کاربری مورد تجزیه‌وتحلیل قرار داده شسده است. نواحی خحشک و نیمه‌عشک در اثر 
بارش cla‏ شدید به‌طور مسلم دارای وقوع شرایط سیلابی شدیدتری هستند. در منطقه موردمطالعه حوضه آبریز شیراز 
با توجه به شرایط اقلیمی منطقه وقوع بارش‌های سیل آسا در زمان‌های کوتاه‌مدت. bo par‏ در ارتفاعات شمالی 
حوضه موحب به راه افتادن حجم عظیمی از دبی گشته است که این مقادیر دبی از طریق مسیلی تحت عنوان رودخانه 
خشک» از Ole‏ شهر شیراز گذر کرده و به سمت دریاچه مهارلو به‌عنوان زهکشی این حوضه سرازیر می‌شود؛ اما با 
روند توسعه شهرنشینی و دست‌اندازی بر حریم مسیل. امروزه شاهد بروز سیلاب و خسارات ناشی از آن هستیم که 
نمونه بارز آن اتفاق فاجعه‌بار سیلاب ۱۳۹۸ در بخش دروازه قرآن شیراز بوده است. لذا ضرورت انجام این تحقیق بر 
مطالعه مقادیر دبی حداکثری بر روند وقوع سیلاب شهری است. اهمیت دوچندان این مطالعه دلالت بر سازوکارها و 


مدریت ble‏ شهری, تغییر کاربری اراضی و عدم رعایت حریم رودخانه بوده است که موجب شده مقادیر حداکثری 


1 Khattak 
2 Kabenge 
3 Zelenakova 


دبی که در زمان‌های گذشته موجب رخداد و خساراتی نشده است. چگونه امروزه با کوچک‌ترین شدت و مدت‌زمان 
کوتاه یک بارش نتواند از طریق مسیل خود زهکشی شده و موجب وقوع سیلاب در حوزه شهری شیراز شود. 
۲- موقعیت منطقه موردمطالعه 

منطقه موردمطالعه حوضه آبریز شیراز دارای محیطی برابر با ۲۳۶/۵ کیلومتر و مساحتی برابر با۱۸۹۸/۲۶ 
کیلومترمربع و محدوده سیاسی استان فارس واقع شده است. این حوضه در مختصات جغرافیایی بین عرض 
جغرافیایی ۲۹ درجه و ۲۰ دقیقه تا ۲۹ درجه و ۵۸ دقیقه شمالی و طول جغرافیایی ۵۲ درجه و ۱۰ adds‏ تا ۵۲ درجه 
OY‏ دقبقه شرقی واقع شده است. حداکثر ارتفاع حوضه از سطح boo‏ ۰ متر و حداقل ارتفاع نیز ۰ متر است. 
خروجی آب از Ob‏ حوضه به سطح اساس دریاچه مهارلو منتهی می‌گردد که تمامی جریانات وارده به این حوضه از 
طریق مسیلی تحت عنوان رودخانه خشک زهکشی می‌گردد و از نظر ژئومورفولوژیکی در واحد زاگرس چین‌خورده 
قرار دارد و از نظر زمین‌شناسی سنگ کف دو ناحیه عمیق در محدوده دشت شیراز که توسط طاقدیس‌های؛ 
سنگ‌آهک‌های آسماری از هم جدا شده است» عمدتاً از رسوبات ماسه‌سنگ. کنگلومرا؛ گچ و Sal‏ به رنگ زرد و 
قرمز تشکیل شده است. رسوبات این حوضه متعلق به دوران چهارم بوده و منطقه وسیعی از دشت شیراز را به‌صورت 
آپرفت پوشانده است. این رسوبات از شن و ماسه و رس تشکیل شده که ضخامت آن حداکثر به ۲۵۰ متر می‌رسد 
(آقا نباتی. ۱۳۸۳). بارش متوسط حوضه طی دوره آماری ۲۰ ساله ۳۱۰,۵۱ میلی‌متر و میانگین دمای سالانه ۱۸,۶۲ 
درجه سلسیوس بوده است و بر مبنای اقلیم آمبرژه بارندگی‌های این ناحیه اکثراً از سیکلون‌های غرب و برخورد هوای 
سرد سیبری و مدیترانه‌ای منشأً می‌گیرد. رطوبت موجود در هوای مدیترانه بر اثر صعود به دامنه‌های زاگرس موجب 
بارندگی در نواحی غربی ایران شده به‌نحوی که میزان آن از غرب به شرق نقصان می‌یابد و شروع فعالیت سیکلون‌های 
مدیترانه‌ای از اوایل پاییز است که تا اواسط بهار ادامه می‌یابد. 


۱۳۲ جغرافیا و مخاطرات محیطی شمارة چهارم 
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شکا \_ نقشه موقعیت حوضه آبریز شیراز e‏ در استان و کشور 


۳- مواد و روش‌ها 

این تحقیق بر OL‏ روش آماری- تحلیلی و از نوع کاربردی است. در این پژوهش از نقشه‌های توپوگرافی 
۰ شیراز-کازرون» زمین‌شناسی ۱/۱۰۰۰۰۰ و داده‌های مرکز مطالعات آب منطقه‌ای استان فارس» در حصوص 
آمار دبی ایستگاه‌های هیدرومتری این سازمان استفاده شده است. در ابتدا محدوده حوضه آبریز شیراز در محیط 
نرم‌افزار Arc Gis‏ معین گردیده است.که حوزه مدنی شیراز نیز در بخشی از این حوضه آبریز واقع شده است. در 
مرحله بعد از داده‌های ایستگاه‌های هیدرومتری oll le‏ مارون‌اقبال آباد. چنارسوخته و پل باغ صفا جهت تعیین 
بیشترین مقادیر دبی این ایستگاه‌ها در یک دوره ٤٤‏ ساله در محیط ۳616 استفاده شده است. پس ازآن از ساختار 
گمبل تیپ یک برای تعیین پیک حداکثر مقادیر دبی ثبت شده در رودخانه شک از طریق رابطه ثابت اویلر- 
ماسکرونی استفاده شده است. همچنین این ساعتار توزیعی از طریق ضوابط جبری. جهت محاسبه گشتاور اول و 
دوم انحراف استاندارد مقادیر دبی» پیش‌بینی‌های جریانات طغیانی تغیبرات چگال از طریق توزیع نرمال گاوس ‏ 
توزیع ارگودیک هولومورفیک حد آستانه» پیش یابی ارگودیک فاصله اطمینان و خروجی‌های مربوطه در محیط 
Graphers 6‏ انجام پذیرفته است و در آخر مقادیر به دست آمده را با استفاده از روش انتگرال جز به جز ریمان» بر 
روی ساختار توزیم گمبل پیاده‌سازی کرده و به بررسی چگونگی روند پیک‌های حداکثری بر وقوع سیلاب ٤٤‏ ساله 


پرداخته شده است. 


1 Gauss 


سال دوازدهم بررسی حداکثر دبی poli‏ سیلاب حوضه آبریز .... ۳ 
در ساختار عمومی توزیع نااریب حد آستانۀ حداکثر گمبل از منفی تا مثبت بی‌نهایت داریم که : 


26- لا‎ 2 Us 


f) = وچ‎ Pec ۰ )( aa 


در توزیع فوق :8 : انحراف استاندارد تعدیل شدۀ گمبل از جامعه از دبی‌های حداکثر و /۸: میانگین تعدیل BLS‏ 
گمبل از دبی‌های حداکثر می‌باشد. همان‌گونه که از تابع گمبل معلوم هست در حقیقت مقادیر انحراف استاندارد و 
میانگین دبی‌های حداکثر حوضۂ شیراز بایاس(اریب) تابع گمبل می‌باشند . جهت رفع نااریبی مذکور ابتدا : ۸ را 
برابر صفر و 7 را برابر یک فرض می‌نماييم . در این صورت جهت محاسبهٌ 6.5.٥‏ در تابع استاندارد گمبل حد 
آستانة حداکثر تیپ یک داریم که : 


x 


رابطه f(X) =e *e“ (Y)‏ 
با توجه به این که تابع احتمال گمبل به‌صورت خودکار چگال می‌باشد؛ Cpl ple‏ داریم که : 


SEE A )۳( رابطه‎ 


t‏ بحث و نتایج 

۱-۶- تحلیل پیک دبی‌های رودخانةٌ خشک شیراز در ساختار گمبل تیپ یک 

یکی از راهکارهای محاسبه در برآورد مقادیر حد آستانه دبی حاصل از بارش‌هاء ساختار گمبل تیپ یک است. 
بدیهی است که با این تحلیل بیشترین و کمترین مقادیر وقوع دبی در حوضة شیراز و Soph‏ وقوع این مقادیر حد 
آستانه هرچند وقت یک بار اتفاق خواهد افتاد را می‌توان محاسبه کرد. در این آنالیزء ٤٤‏ سال دیتاهای حداکثر دبی 
رودخانهةٌ خشک مورد استفاده قرار داده شده است. جهت رفع بایاس از مقادیر حداکثری دبی حاصل از بارش‌های ££ 
ساله‌ی حوضة شیراز ایستگاه‌های علی‌آباد مارون اقبال‌آباده چنار سوخته و پل باغ صفا واقع در شهر شیران برداشت 
گردیده است. از داده‌های هیدرومتری ایستگاه پل فسا که در محل اوت‌لت (خارجی) رودخانه خشک به این 
رودخانه می‌پیونده صرف‌نظر شده است. به دلیل این که نقشی در توزیع روان آب احتمالی بر سطح حوضه و شهر 
شیراز ندارد. بدیهی است درصورتی که در مطالعات کل حوضۀ هیدروژئومورفولوژیک مهارلو» مورد بررسی قرار گیرد 


نه‌تنها این شاخه از رودخانه, بلکه شاخه‌های جنوبی حوضه شامل زیر حوضه سروستان بوده و Ub‏ در محاسبات 


۱۳ جغرافیا و مخاطرات محیطی ls‏ چهارم 


طریق at‏ چگال احتمال رابطه )£( داریم که : 


F(X) = (f? ()(e e “dx)) = 0.577215664902 )٤( رابطه‎ 


این تابع احتمال برابر است با 0.577215664902 که این مقدار برابر با ثابت اویلر - ماسکرونی بوده و حاصل 


مذکور از محاسبه انتگرال جز به جز از ساختار ریمان در نمودار زیر قابل‌محاسبه می‌باشد. 


Euler—-Mascheroni constant 


شکل -Y‏ مقدار گشتاور اول توزیع حد آستانة نوع اول گمبل 
در ادامه جهت رفع بایاس » میانگین و انحراف استاندارد دبی‌های رودخانة خشک در تابع احتمال نوع یک گمبل 
با محاسبه‌ی گشتاور دوم توزیع حد Seal‏ گمبل از ساختار جز به جز ریمانی رابطه (۵) داریم که : 


+ 


F(X^2) = (J S (x2)(e*e-© “dx)) = 1.97811199066 (0) رابطه‎ 


این مقدار فوق» سطح محصور زیر منحنی انتگرال جز به جز گشتاور دوم گمبل» از ساختار ریمان را نشان داده 
است. نمودار زیر از طریق رفع بایاس مقادیر حداکثر سالانه دبی Seog)‏ خشک حوضه مطالعاتی به دست آمده 


ual 


سال دوازدهم بررسی حداکثر دبی مقادیر سیلاب حوضه آبریز ... ywo‏ 


D 
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Variance of Gumbel Distribution 
شکل ۳- مقدار گشتاور دوم گمبل از ساختار ریمان‎ 


VAR(P) = 1.97811199066 — 0.577215664902? = (A) رابطه‎ 
STD(P) = (V) رابطه‎ 4 


4 (1.97811199066) — (0.577215664902?) = 1.28254983016 


با محاسبة گشتاور اول و دوم و انحراف استاندارد تابع احتمال چگال گمبل» در رفع بایاس دبی‌های حداکثری 
حوضة رودخانةٌ خشک شیراز از طریق محاسبة ضریب کشیدگی گمبل تیپ یک روابط (۸ و (A‏ به دست آمده‌اند. 


رابطه (Sn) = 0.450053207457 W)‏ را 


Geometric Gumbel Kurtosis( 
STD(P) 1.28254983016 


Gumbel Ratio (STDP) = (a) = (zzz) = 0.779696801235 = ((Bm)/V(6)) (4) رابطه‎ 


STD(P) 1.28254983016 


Mean = (J? (x)(e*e~® “dx)) = 0.577215664902 (\+) رابطه‎ 


با این شرط که در حالت استاندارد ۸۸ برابر صفر و / برابر یک فرض می‌گردد. مقدار Mean‏ برابر خواهد بود با 
cal‏ اویلر -ماسکرونی » یعنی مقدار 0.577215664902. به همین طریق ply MSI‏ صفر گردد بنابراین Mode‏ 
در ساختار توزیع یک گمبل برابر با رابطه (۱۱) : 


۱۳۹ جغرافیا و مخاطرات محیطی els‏ چهارم 


Mode =u =0 )۱۱( رابطه‎ 


و درصورتی که در حالت استاندارد گمبل تیپ یک .مقدار ۸4 برابر phe‏ فرض شود در محاسبه کمیت مقدار مد 
توزیع برابر 0,366512920582 است, 


fX) = (e%e-8*)~ (eee) = >-(0)-و)‎ ۳ ( = 0.366512920582 (Y) رابطه‎ 


O= EE a) S05 )۱۳( رابطه‎ 
F(X) = Median = u — B(In(In(2)) <- 9 )۱۶( abl, 


راست چولگی دارد که با توجه به نوع عملکرد این ساختار کاملاً منطقی است. هم‌چنین میزان کشیدگی جبری توزیع 


نیز ۵,۶ می‌باشد . 


).368 Mode Gumbel =0.36787944117144. 3 
Median Gumbel 20.6 


Mean Gumbel = 0.3202 


).276 | 


(.184 7 
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Standard Distribution of Gumbel Maximum Extreme Value Type 1 


شکل ۶- میزان کشیدگی جبری توزیع گمبل تیپ یک 
در محاسبة توزیع تجمعی احتمال حد آستانۀ گمبل تیپ یک رابطه (۱۵) برقرار است که : 


P = ee” )۱۵( رابطه‎ 


سال دوازدهم بررسی حداکثر دبی مقادیر سیلاب حوضه آبریز ... aii‏ 


در شکل ۵ محور افقی سری زمانی و محور عمودی احتمال وقوع مطلق است. Lb‏ در نظر داشت که محاسبۀ 
فوق و نمودار زیر هردو در حالت کلاسیک توزیع تجمعی را نشان می‌دهند که در تحلیل‌های حوضوی Lb‏ بایاس 
در دیتاهای بەدستآمده رفع 
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Cumulative Distribution of Gumbel Maximum Extreme Value Type 1 


شکل ۵ - حالت کلاسیک توزیع تجمعی سری زمانی و احتمال وقوع مطلق 


قابل محاسبه است : 


اد 


(- 


در واقع ساختار فوق که در تحلیل آن نااریپی‌های 8 و / اعمال می‌گردد نه‌تنها باید در محاسبات دبی‌های حد 
آستانة حداکثر تیپ یک گمبل حوضء شیراز Gb‏ (جدول ۱) اعمال شود بلکه در pls‏ آنالیزهای مشابه نیز Gale‏ 


است. 


۱۳/۸ جغرافیا و مخاطرات محیطی ishi‏ چهارم 


جدول ۱- مقادیر دبی 26 ساله حوضه هیدروژئومورفیک شیراز 


CV =STDP/MEAN G.D.R (G.D.R)-CV‏ انحراف استاندارد جامعه ي متحرک میانگین متحرک سورت صعودی حداکثر دبی روزانه در سال به متر مکعب سال 
0.497038512 0.39894228 0.895980792 38.48801561 42.95629545 0.248 1387 
2 0.39894228 0.873081302 38.37149686 43.94951163 0.669 1386 
8 0.39894228 0.850001158 38.23305208 44.98 1.23 1388 
0.427897688 0.39894228 0.826839968 38.07356052 46.04707317 24 1355 
0.405245991 0.39894228 0.804188271 37.90802777 47.13825 5.77 1360 
4 0.39894228 0.784254474 37.80026128 48.19897436 8.18 1357 
5 0.39894228 0.765965555 37.72541613 49.25210526 9.67 1392 
0.349415386 0.39894228 0.748357666 37.65877356 50.32189189 11.2 1358 
0.332484186 0.39894228 0.731426466 37.60161877 51.40861111 12.79 1375 
0.316245718 0.39894228 0.715187998 ,37.55595215 52.512 12.9 1396 
4 0.39894228 0.698169174 37.47566783 53.67705882 15.01 1359 
0.283308584 0.39894228 0.682250864 37.42063292 54.84878788 15.32 1354 
7 0.39894228 0.665650867 37.33240481 56.0840625 16.6 1398 
1 0.39894228 0.649243701 37.23915263 57.35774194 17.55 1353 
0.233706278 0.39894228 0.632648558 37.12676974 58.68466667 19.8 1389 
0.218135357 0.39894228 0.617077637 37.04040437 60.02551724 21.09 1361 
0.202606536 0.39894228 0.601548816 36.94476508 61.41607143 24.6 1379 
0.189215928 0.39894228 0.588158208 36.92435445 62.77962963 25.8 1378 
0.175726212 0.39894228 0.574668492 36.8948234 64.20192308 27.5 1393 
0.162551611 0.39894228 0.561493891 36.87330384 65.67 28.57 1363 
0.149014959 0.39894228 0.547957239 36.83140247 67.21583333 28.63 1370 
0.133980647 0.39894228 0.532922927 36.71491411 68.89347826 28.9 1385 
0.117434103 0.39894228 0.516376383 36.51367817 70.71136364 31.7 1372 
0.101863768 0.39894228 0.500806048 36.34301792 72.56904762 35 1394 
0.087735812 0.39894228 0.486678092 36.2396722 74.4475 35.4 1391 
0.071879607 0.39894228 0.470821887 36.0191332 76.50263158 35.8 1390 
0.053921632 0.39894228 0.452863912 35.6693282 78.76388889 36.88 1369 
6 0.39894228 0.433147076 35.18351779 81.22764706 38 1364 
0.012274504 0.39894228 0.411216784 34.51316768 83.929375 48.5 1384 
0.0006- 0.3989 0.39830371 34.37015815 86.291333 51.84 1368 
0.01274168- 0.39894228 0.3862006 34.27613084 88.75214286 53.7 1381 
0.02595292- 0.39894228 0.372989363 34.10930346 91.44846154 57.6 1362 
0.03812235- 0.39894228 0.360819928 34.01419396 94.26916667 58.81 1367 
0.05300804- 0.39894228 0.345934237 33.72607222 97.49272727 61.38 1374 
0.0697495- 0.39894228 0.329192776 33.28270638 101.104 70.5 1395 
۰-۰6 ۰ 0.39894228 0.31955068 33.39446631 104.5044444 75.1 1371 
0.08779334- 0.39894228 0.311148938 33.6600921 108.18 82 1382 
0.09163204- 0.39894228 0.307310242 34.3941623 111.92 83.94 1373 
0.09837375- 0.39894228 0.300568531 35.04128121 116.5833333 91.7 1377 
0.09951225- 0.39894228 0.299430028 36.39871426 121.56 93 1376 
0.1080922- 0.39894228 0.290850084 37.43240575 128.7 93.7 1380 
0.13973478- 0.39894228 0.259207502 36.38409298 140.3666667 94.1 1366 
0.27967622- 0.39894228 0.119266055 19.5 163.5 144 1383 
۰-8 0.39894228 0 0 183 183 1397 


بنابراین می‌توان مشاهده نمود که پیک دبی‌های ably‏ شده طی ٤٤‏ سال از ۰,۲۶۸ مترمکعب بر انيه تا ۱۸۳ 
مترمکعب در ثانیه در سال ۱۳۹۷ متغیر بوده است. بدیهی است که حتی دبی‌های پیک حداکثری در سال‌های 
آماربرداری نیز نمی‌توانند Ge‏ روان آبی بیش در نظر گرفته شوند که سبب ایجاد = ob‏ طغیانی در رودخانه خشک 
شده است؛ اما امروزه بنا به دلیل مدیریت‌های غلط شهری و عدم رعایت حریم رودخانه. تبدیل به سیلاب شهری 
گشته است. در تعیین دبی که می‌تواند به‌عنوان Liv‏ جریان طغیانی در رودخانةٌ خشک در پیش‌بینی آینده شود از 
طریق اعمال ضریب تغییرات چگال توزیع نرمال گاوس قابل محاسبه است که بر اساس روابط (۱۷ و OA‏ برابر است 
L‏ 


10) = SS 2 dx = Ven )۱۷( رابطه‎ 


سال دوازدهم بررسی حداکثر دبی مقادیر سیلاب حوضه آبریز ... ra‏ 


1 


1 


Gauss Density Ratio = 


پس باید ضریب تغییرات دبی. حداقل مقدار 0.39894228 باشد که همان مقدار چگال توزیع گاوس است؛ 
بنابراین بر اساس محاسبات تحلیلی فوق‌الذکر از دبی‌های ٤٤‏ ساله قابل‌تحلیل در ساختار گمبلء آن‌هایی روایی دارند 
که مقدار آن‌ها حداقل بیشتر از 51.48 مترمکعب در acl‏ باشد. 

در ادامه برای پیش‌بینی فاصلة اطمینان دپی حداکثر در محاسبۀ حد آستانه دبی‌های رودخانه (SLES‏ ابتدا باید 
بایاس‌های ساختار گمبل محاسبه گردد . بنابراین در یک ساختار مجموع کلاسیک نامعین در برآورد فاصلٌ اطمینان 
بین دو مقدار از دبی‌های حد آستانة رودخانه خشک شیراز برابر توزیع حد آستانه‌ای حداکثر گمبل تیپ یک رابطه 
(VY)‏ حاکم است: 


رابطه (۱۹) اق ق F(X) = (f‏ 


و چون al‏ هولومورف فوق به‌صورت هندسی چگال می‌باشد بنابراین کافی است که مقادیر ( 0 , 4) به توزیع 
احتمال مذکور تزریق گردد. نباید از نظر دور داشت که توزیع احتمال گمبل تیپ یک از یک شبکه ارگودیک پیروی 
می‌نماید. در این شبکه مقادیر × از لگاریتم طبیعی روی محور عرضی (Y)‏ محاسبه می‌گردد و سپس این سیستم 
ارگودیک. به‌صورت نمای ثابت )©( قرار می‌گیرد (شکل CV‏ بنابراین می‌توان به مکانیزم عملگر متغیرهای سه‌گانة 


(a,b,x)‏ در این شبکة ارگودیک هولومورف به‌حوبی پی برد. 


The critical region in the Exterm Gumbel Type 1 
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شکل 1 ساختار حد Gul‏ استاندارد تیپ یک گمبل 


۱۳۰ جغرافیا و مخاطرات محیطی els‏ چهارم 


در توزیع استاندارد گمبل اریب (Bias)‏ به‌راحتی ابتدا مقادیر (a,b)‏ محاسبه می‌گردد . این دو مقدار یک سری 
زمانی هستند که تابع سیر لگاریتم طبیعی هستند. به بهترین وجهی ساختار ارگودیک گمبل با رفع بایاس این دو مقدار 
را بهضورت aly fae a asl‏ کر ان hots‏ ارگردیک گم اھا Al pal‏ انتاندارد یک دی هی 
روزانة رودخانة حشک را در ضریب 0.779696801235 ( تابع شسمارة 4( مقیاس می‌نماید. شکل ۷ دارای 
E ao Sioa Glass‏ ا ھ See‏ وک ی all‏ قاد See Sl sabes‏ 
دبی‌های حداکثر روزانه رودخانة خشک. اگر به‌صورت مستقیم ( روی متغیرهای 0ب۵, ) یا به‌طور غیرمستقیم فقط 
روی متغیرهای اه اعمال گردد. بایاس انحراف استاندارد گمبل در این برآورد رفع خواهد شد. همچنین توزیع حد 
آستانة گمبل chlo‏ یک kurtosis‏ هندسی است که مقدار آن برابر با 0.450053207457 است. dlals)‏ شمارة 


AA 


QA3&P%&R in Year 


4 3 -2 1 0 1 2 3 4 5 6 7 8 9 10 Y 
N n n q q q q q 
n Q اب‎ bed = ها‎ Q bad © 
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Return Period (Year) 


شکل ۷- انحراف استاندارد دبی‌های حداکثر روزانه 


از یک تابع‌نمایی مضاعف هولومورف هستند باید در alag‏ اول با مقدار تابع شمارۀ ٩‏ و سپس با مقدار al‏ شمارة ۸ 


FO) = fee OP) yay (1+) رابطه‎ 


سال دوازدهم بررسی حداکثر دبی مقادیر سیلاب حوضه آبریز ... ۳۱ 


باید در نظر داشت که مقدار تابع دیفرانسیل شمارهةٌ 20 بر خلاف تابع احتمال شمار( ۳) به‌صورت خودکار چگال 
os ۳ A ۳ ۳ Z ert 4 ۳ ۲ :‏ 
نمی‌باشد. همچنین تابع باید در مقدار (GB)‏ ضرب گردد تا به کورتیس مورد نیاز چگال برسد. تفسیر عملی فوق روی 
شکل ۷ توزیع حد آستانة ماکزیمم تیپ یک گمبل می‌باشد. ورای اعمال دوره بازگشت دبی‌های حد آستانة بین ۱6۰ 


150 = 86.29133 + )34:37016))0.779697 *۷( — 0.450053208) (1) abel, 
250 = 86.29133 + (34.37016((0.779697 * ¥) — 0.450053208) )۷۲( bal 
¥(150m3) = (Snn + 0.450053208 |) /0.77969 = 2.95455 
MISU 34.37016 شش‎ me 
Y (250m?) = (= 5 ات‎ + 0.450053208 |) /0.77969 = 6.68619 
eS 34.37016 0 


در دو dls‏ ۱ و ۲۲ ضرایب بایاس از ضوابط A‏ و٩‏ استخراج شده است.پس با محاسبة مقادیر Y‏ در توزیع 


FOE € = 0.04952 T 


در تابع اولیه احتمال فوق فاصله اطمینان احتمال وقوع دبی‌های بین ۱۵۰ تا ۲۵۰ متر کعب در انیه حداکثر برابر با 
٩‏ درصد است و باید کاملاً دقت داشت که این کمیت یک Lob‏ اطمینان روی توزیع استاندارد تابع احتمال گمبل 
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شکل ۸- فاصله اطمینان مقادیر توزیع استاندارد تابع احتمال گمبل تیپ یک 


۱۳۲ جغرافیا و مخاطرات محیطی els‏ چهارم 


در حقیقت مقادیر Y‏ های به‌دست‌آمده روی مقادیر مورد پیشنهاد شده از نمودار A‏ به نواحی روی نمودار ۷ 
انتقال یافته‌اند . در این تحلیل حتماً و حتماً باید مقادیر Y‏ بین 6- و © واقع گردد. در غیر این صورت فرض 11 به 


وقوع خواهد پیوست که رد فرض HO‏ می‌باشد؛ به‌عبارتی‌دیگر جهت قبول فرض 10 باید رابطه (VE)‏ حاکم باشد : 


۲۱-9 ۷ > (YE) رابطه‎ 


مقادیر ۲ ساختار ارگودیک لگاریتم نپرین» منشعب از یک سری زمانی معکوس است. به این معنی که در سری 
زمانی گمبل عامل عملگر به جای ‏ که Sl‏ احتمال وقوع واقعی است. از 4 استفاده شده است. 


رابطه (۲۵) << 


پاک وه کشت کزان به سال ی E‏ بت این کم bagels ue‏ کر فک نله بانه شاف یک 
شال و کسری از یک سال باشد. در محاسبات حد آستانه ماکزیمم تیپ یک گمبل؛ از احتمال عدم وقوع استفاده 


q=1-p (VV) رابطه‎ 


پس در محاسبةٌ ۲ داریم که : 


Y = —In(—In(q)) (YV) رابطه‎ 
بنابراین:‎ 


Y =—In(—In (1 ۳ Ə) (YA) رابطه‎ 


به خوبی آشکار است که عامل ۷ در ساختار تیپ یک گمبل» یک عملگر ارگودیک از یک سری زمانی است 
که ترکیبی از احتمال وقوع و عدم احتمال وقوع (08-0)هست. بنابراین جهت محاسبه مقدار دبی حداکثر با استفاده از 
مقدار Y‏ محاسبه شده در ساختار توزیع تیپ یک ماکزیمم برابر است با : 


Q/M*3 = 0 + (a(Q) = ((0.779696801235 *Y)— 0.450053207457)) )۲۹( رابطه‎ 


بر این اساس علت ارگودیک بودن ساختار گمبل تیپ یک ماکزیمم» این است که در این سیستم پژوهشگر با 
حداقل ده بعد شامل احتمال وقوع عدم احتمال وقوع, محور طولی که با ۷ ارگودیک می‌گردد. محور ۷ که با 
P‏ مطلق و نسبی ارگودیک می‌گردد و نموداری که یکی دارای ضریب هندسی Skewness‏ برابر 0.779 مثبت به 
راست بوده» با kurtosis‏ در حد 0.45 هندسی است. 

در تحلیل بر روی ob‏ اطمینان بین دو مقدار ارگودیگ دوره بازگشت. احتمال و مقدار دبی حداکثر بیانگرء 
ساختار حداکثرهای آستانه در هر نقطه بوده و این ساختار کاملا ار گودیک است. که این ساختار با ساختارهای 
کلاسیک احتمال دامنه به راست» هم‌خوانی داشته» وضمنا عمل کرد آن بسیار مجزاست. اعمال این ساختار روی پیک 
دبی‌های روزانه رودخانه حشک شیراز در طول یک سال طی دوره آماری £6 ساله در شکل ٩‏ قابل تفسیر ومشاهده 


است. 


433,34 95:14 98-19 99.33 99.75 99.91 99.97 99.98 99.99 ۲ 


244 151.2 178 204.8 231.6 258.4 285.2 312 338.8 ۶ 


20.6 55.1 148.9 403.9 1097.1 2981.5 8103.6 22027 T/yeat 


Po o o o oo 659 ag 79 


Shiraz Rodkhaneh KhoskCumlative Distribution of the Gumbel maximum Extreme Value 


شکل -٩‏ ساختار حد آستانه تجمعی حداکثر تیپ یک گمبل در پیش بینی فاصله اطمینان دبی 


که در تفسیر احتمال کلاسیک راست دامنه داریم که : 
پذیرش فرض 90 > 9 :۱10 
رد فرض 60 > :110 
در بیان احتمال کلاسیک فوق برای مثال در رودخانه خشک شیراز . دبی پیک ۱۲۶ مترمکعب روزانه طی یک 
ellis‏ اک بتر از ٤‏ درصد. طی ۸ سال احتمال وقوع ندارد و این دقیقاً از ساختارهای نرمال و گاما بیس تبعیت 


می‌نماید. Lb‏ توجه نمود که در تمام تحلیل‌های گمبل تیپ یک اگر ضریب انحراف استاندارد و چولگی گمبل اعمال 


۱۳ 


جغرافیا و مخاطرات محیطی 


els‏ چهارم 


گردد نتایج نااریب Un Bias‏ و در غیر این صورت اریب یا Bias‏ خواهند بود. تبیین محاسباتی بر روی دبی‌های 


روزانه حداکثری ably‏ شده در رودخانه خشک شیراز که شکل ٩‏ از آن استخراج شده است به شرح جدول شمارة 


است. 


T/Year 


رط طط 


1 
1.000005 
1.000618 
1.011444 

1.07065 
1.238077 
1.581977 
2.198958 

3.24887 
5.000268 
7.900331 
12.68933 
20.58969 
33.61797 
55.09968 
90.51806 
148.9137 
245.1923 

403.929 
665.6418 
1097.133 
1808.542 
2981.458 
4915.269 
8103.584 
13360.23 
22026.97 


۵- نتیجه گیری 


جدول ۲- تبیین محاسباتی دبی‌های روزانه حداکثری رودخانه خشک 


96۳ 
3.50739E-63 
8.05483E-38 
1.94234E-22 

4.1509E-13 
1.89218E-07 
0.000511929 
0.061797899 
1.131428638 
6.598803585 
19.22956455 
36.78794412 
54.52392119 
69.22006276 
80.0010713 
87.34230185 
92.11936552 
95.14319929 
97.02540026 
98.18510731 
98.89524805 
99.32847021 
99.5921568 
99.75243174 
99.84976904 
99.90885337 
99.94470686 
99.96645936 
99.97965523 
99.98765978 
99.9925151 
99.99546011 


P 
3.50739E-65 
8.05483E-40 
1.94234E-24 

4.1509E-15 
1.89218E-09 
5.11929E-06 
0.000617979 
0.011314286 
0.065988036 
0.192295646 
0.367879441 
0.545239212 
0.692200628 
0.800010713 
0.873423018 
0.921193655 
0.951431993 
0.970254003 
0.981851073 
0.988952481 
0.993284702 
0.995921568 
0.997524317 
0.99849769 
0.999088534 
0.999447069 
0.999664594 
0.999796552 
0.999876598 
0.999925151 
0.999954601 


Q/M/*3 


3.827177 
17.22633 
30.62548 
44.02463 
57.42378 
70.82293 
84.22208 
97.62124 
111.0204 
124.4195 
137.8187 
151.2178 
164.617 
178.0161 
191.4153 
204.8144 
218.2136 
231.6127 
245.0119 
258.4111 
271.8102 
285.2094 
298.6085 
312.0077 
325.4068 
338.806 


STD(p) 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 
34.37016 


AVE 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 
86.29133 


G.G.K 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 
0.450053 


B 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 
0.779697 


۲ 


موجب وقوع این اتفاق می‌شود. ازجمله شرایط توپوگرافیک تغییرات اقلیمی. ضعف مدیریت‌های حوضوی جهت 


اندازه‌گیری مقادیر دبی حداکثری ایستگاه‌های هیدرومتری واقع شده بر روی رودخانه خشک بر ایجاد سیلاب در 


حوضه آبریز شیراز طریق محاسبه ضوابط جبری توزیع گمبل تیپ یک برای بررسی روند یک دوره ٤٤‏ ساله پرداخته 


شد. بدین‌صورت که ابتدا جهت رفع بایاس میانگین و انحراف استاندارد مقادیر دبی و محاسبه گشتاور اول و دوم از 
ساختار انتگرال جز به جز ریمانی» استفاده شده است. تا توزیع استاندارد حد استانه ماکزیمم گمبل تیپ یک» در 


تحلیل دبی‌های حداکثری به دست آید که برابر با VOW OTEA‏ بوده است. شرط لازم مقدار میانگین برابر با 
ثابت اویلر- ماس‌کرونی بوده و نشان‌دهنده نوع عملکرد منطقی توزیم جبری است. در مرحله بعد با رفع بایلاس 
ھان aa‏ یت EE‏ که OE al: a‏ ا WAP WAN a‏ مکی ر ا کی ردو که ای 
مقادیر امروزه تحت عنوان سیلاب شهری مشکل‌ساز شده درصورتی که در بستر اصلی حریم رودخانه هرگز سیلابی» 
بهحساب نیامده و نهایتاً یک = Ob‏ طغیانی بوده است. همچنین در پیش‌یابی فاصله اطمینان دبی حداکثری برای 
اعمال دوره‌های بازگشت در ساختار ارگودیک هولومورفیک گمبل» مقادیر دوره بازگشت دبی‌های حد آستانه بین 
۰ تا ۲۵۰ مترمکعب بر انیه در طول سال» برابر با ٤٩‏ درصد طبق (حدول شماره ۲ ) بوده است. به‌طوری‌که به 
وقوع پیوستن دبی پیک طی صدسال با احتمال ۷۰ در صد با حدود ۲۰۰ مترمکعب در ثانیه امروزه موجب سیلاب 


گردیده و اگر این بستر به تصرف مهندسی غلط شهری در نیامده بود به معنای وقوع سیلاب تلقی نمی‌گردید. 
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